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In Bose supefluids, the single-particle Green's function can be directly related to the superfiuid 
velocity-velocity correlation function in the hydrodynamic regime. An explicit expression for the 
single-particle spectral density was originally written down by Hohenberg and Martin in 1965, 
starting from the two-fluid equations for a superfiuid. We give a simple derivation of their results. 
Using these results, we calculate the relative weights of first and second sound modes in the single- 
particle spectral density as a function of temperature in a uniform Bose gas. We show that the 
second sound mode makes a dominant contribution to the single-particle spectrum in relatively high 
temperature region. We also discuss the possibility of experimental observation of the second sound 
mode in a Bose gas by photoemission spectroscopy. 

I. INTRODUCTION 

In the low-frequency dynamics of superfiuid Bose and Fermi atomic gases, the most dramatic effects related to 
superfluidity are described by Landau's two-fluid hydrodynamics analogous to the case of liquid 4 He [lj]. These 
equations only describe the dynamics when collisions are sufficiently strong to produce a state of local thermodynamic 
equilibrium [2|. This requirement is usually summarized as lot <C 1, where w is the frequency of a collective mode 
and r is the appropriate relaxation rate. In this regime two sound modes can be distinguished: the first sound 
mode consists of an in-phasc oscillation of the superfiuid and normal fluid components, while the second sound mode 
consists of an out-of-phase oscillation of the superfiuid and normal fluid components. The occurrence of two distinct 
modes arises from the presence of both a superfiuid component and normal fluid component, which are coupled to 
each other. Recently, there has been renewed interest in second sound mode in superfiuid Bose and Fermi gases ^-Q|- 
The study of ultracold gases in collisional hydrodynamic regime has been difficult because the density and the s-wave 
\ scattering length are typically not large enough. However, Feshbach resonances allow ones to achieve conditions where 
the Landau two-fluid description is correct. Recent experiments have begun to observe sound propagation in trapped 
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superfiuid Fermi gases with a Feshbach resonance 



At unitarity, the magnitude of the s-wave scattering 



length that characterizes the interactions between fermions in different hyperfinc states diverges (\a s \ — > oo). Owing 
to the strong interaction close to unitarity, the dynamics of superfiuid Fermi gases with a Feshbach resonance at finite 
temperatures are expected to be described by Landau's two-fluid hydrodynamic equations HQ. 

More recently, sound propagation in a Bose-condensed gas in a highly elongated (cigar-shaped) trap has been 
observed in Ref. , where the thermal cloud is in the hydrodynamic regime and thus the system is described by the 
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two-fluid model. In this experiment, the sound wave in a highly elongated trapped gas can be excited by a sudden 
modification of a trapping potential using the focused laser beam. The resulting density perturbations propagate 
with a speed of sound. This experimental work has reported some success, with evidence for a second sound mode in 
superfluid Bose gases, but first sound mode was not observed. In order to clearly demonstrate the two fluid dynamics 
of superfluids, it will be important to observe both first and second sound modes. In the case of a sound wave excited 
by a sudden modification of a trapping potential, the thermal density perturbations (first sound) is so small that one 
cannot distinguish small density perturbations from signal-to- noise in the thermal cloud in this regime 

In general, the single-particle Green's function of Bose superfluids in the low-frequency, long-wavelength regime 
is directly related to the superfluid velocity- velocity correlation function ll| , which is therefore coupled to first and 



second sound modes. Making use of this exact relation, In this paper, we propose that first and second sound 
can be probed by measuring the single-particle spectral density. This type of quantity is directory related to the 
tunneling current spectroscopy jl^ . Resent experience on ultracold Fermi gases by JILA group [13J showed that 
the momentum-resolved photoemission-type spectroscopy is a powerful technique to directly probe single-particle 
excitations of ultracold atomic gases. 

In Sec. II, we discuss the relation between superfluid velocity-velocity correlation function and the single-particle 
Green's function of Bose superfluids, which was discussed by Hohenberg and Martin (which will be referred to as 
"HM") llj]. We also give a simple derivation of the explicit expression for the single-particle Green's function in the 
two-fluid hydrodynamic regime, following the approach analogous to the derivation of the density correlation function 
given in Ch.14 of Ref. fl. 

In Sec. Ill, we calculate the relative weights of first and second sound mode in the single-particle spectral density. 
For this purpose, we use the Hartree-Fock-Bogoliubov (HFB)-Popov approximation [ijj], for calculating various ther- 
modynamic variables. We also compare the relative weights of first and second sound modes in both the single-particle 
spectral density and the dynamic structure factor 

In Sec. IV, We show the single-particle spectral density in connection with the rf-tunneling current spectroscopy 
12 , [ljj [rH . We will show that both first and second sound mode can be observed by photoemission-type spectroscopy. 



For comparison, we also show the he single-particle spectral density in the collisionless limit using the HFB-Popov 
approximation. 

II. SUPERFLUID VELOCITY- VELOCITY CORRELATION FUNCTION AND THE 
SINGLE-PARTICLE GREEN'S FUNCTION 

According to HM theory, the superfluid velocity-velocity correlation function in the non-dissipative hydrodynamic 
limit is given by 11 1 
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where u\ and U2 are first and second sound velocities, which satisfy u\u\ 
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and normal fluid densities and s is the entropy per unit mass. The explicit expression for A„ s Vs is given by 



Here p s and p n are superfluid 



dP\ T s ( dP\ Ts 2 
2 



dp J s c v p \dT J c 



(2) 



3 



In the low-frequency region, the single-particle Green's function G(q, to) is related to the superfluid velocity- velocity 
correlation function through 11 1 

,2 



G(q,w) = — — x« 3 ,^(q>w), 



(3) 



where no is the condensate density. The single-particle spectral density A(q,uj) is related to the Green's function 
G(q, ui) through 

dui A(q.uj) 



This gives 
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where X\ and X 2 are defined by 
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HM gives a systematic way to calculate various correlations functions in uniform superfluids in the two-fluid hydro- 
dynamic regime. However, their detailed derivation is quite involved, which closely follows the the earlier paper by 
Kadanoff and Martin 17[ on a normal fluid. In this section, we give an alternative derivation, analogous to Sec. 14.3 



of Ref. [l| for the density-density correlation function. We start with the non-dissipative two-fluid hydrodynamic 
equations: 
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The total mass density and mass current are given by the sum of two components 



mn = p = p s + p r , 



(10) 
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mj = p s v a + p„v n . 



(12) 



To calculate the superfluid velocity- velocity correlation function, we add a time-dependent external current 5j ex (r, t) 
that is only coupled to the superfluid velocity v s . That is, the continuity equation becomes 
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where as the entropy equation (jTUJ) is unchanged. 
Taking time derivative of (fl~3"l) Cwc obtain 
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Taking time derivative of (|10[) and linearize it in fluctuations, we obtain 
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This is the same as in the case without the external current, since the derivation does not involve the continuity 
equation. It can be rewritten in terms of the local entropy per unit mass s = s/p. Using 



we see that 



Using dHJ) and (32) in (JT9J), we obtain 
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Let us now use P and T as independent variables. We thus express fluctuations of p and s in terms of P and T as 
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We then obtain 
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We consider an external current that excites modes of frequency uj and wavevector q, namely 
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and look at the plane- wave solutions of Eqs. (|23l) . (l24|) 

SP(r, t) = SP^e^-^ , ST(r, t) = ST^e^-^ . 
This gives two coupled algebraic equations 
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We note that the determinant of the coefficient matrix is given by 
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In obtaining the final expression, we have made used of the following relations: 
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One can also write the determinant in the compact form 
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We have thus obtained the expressions for the temperature and pressure fluctuations. Using the Gibbs-Duhem relat' 

nd/j, = 5P — nsdT, (3 
we can derive the expression for the fluctuation of the chemical potential. 
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In deriving the final expression, we have used Eqs. f|29[) - (f3T|) and analogous formulas. The superfluid velocity can be 
then obtained from the linearized form of (j9|), i.e. 



imujSvsq,^ = —iqSfi 



(36) 



On the other hand, the superfluid velocity can be written in terms of the superfluid velocity-velocity correlation 
function as 



£v sq , w = mxi, s ,^(q, w) — q • <5j cx (q, w). 



Using the relation (|36|) . one can write ([37]) as 
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q 2 

Comparing (|38[) with (|35p . we obtain the expression (JT|) for Xv 3 .v a , where X Vg . Vs is now given by 
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Finally, using the thermodynamic identity 
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we can show that ([39)) agrees with the HM expression ([2]). 



III. AMPLITUDE OF FIRST AND SECOND SOUND 



A. HFB-Popov approximation 
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In order to calculate the single-particle Green's function explicitly, we must specify a microscopic approximation for 
calculating thermodynamic various variables. Here we use the Hartrce-Fock-Bogoliubov (HFB)-Popov approximation 
[14) 1 . For a uniform Bosc gas, HFB-Popov approximation gives the quasiparticle excitation spectrum 
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Here, e® = is the single-particle energy of a noninteracting gas. The noncondensate atom fraction is given by 
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Here, f(E) = \ j{eP E — 1) is the Bosc distribution function. The quasiparticle amplitudes u q and v q are given by 
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Solving (|4Tj) and (|42j) self-consistently, we obtain the condensate atom number Nq, noncondensate atom number N, 
and quasiparticle energy spectrum E q . The dimensionless interaction parameter is defined by 

gn 
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FIG. 1: Quasiparticle energy spectrum for g' = 0.5 at the temperature T = 0.5T°.The solid line is the quasiparticle excitation 
spectrum Eq in Eq. (|41|) and the dashed line is the sound-like energy spectrum E q = chq with c = */ 1 gna/m 



were T® is the BEC transition temperature of an ideal Bose gas. In Fig. [TJ we plot the quasiparticle energy spectrum 
for the interaction g' = 0.5 at the temperature T = 0.5T®. Here the wavenumber is normalized in terms of the healing 
length £ = -yg== ■ We will calculate the thermodynamic functions using the results obtained above. 
In the HFB-Popov approximation, the thermodynamic potential is given by 

fi = —fxrioV + \gn\V + fe B T^ In (l - e"^) + £ (45) 

The pressure is the given by P = — f2/V. The entropy is given in terms of the quasiparticle excitation spectrum as 

S = k B J2{[l + f(E q )]\n[l + f(E q )]-f(E q )\nf(E q )}. (46) 

Using these thermodynamic quantities, one can calculate the sound velocities. The first and second sound velocities 
are given by 

= ± \Z(^r^) 2 -^ (47) 

where C? = (f C$ = (g)^ C| = Cj - C| = (|^)^. 

Using the above results, we calculate the relative weights of first and second sound in A(q, uf) (see Eq. [6]). We take 
the interaction parameter from the experiment of Rcf. which reports the observation of second sound. In this 
experiment, total number of 23 Na atoms N = 1.7 x 10 8 , radial trap frequency w Ia d/2w — 95Hz, and the aspect ratio 
w ra d/wax ~ 65. We estimate the interaction parameter for a uniform gas using the average density of the trapped gas, 
and obtain g' ~ 0.5 (na 3 ~ 0.07). In Appendix B, we evaluate the characteristic collisional relaxation time and confirm 
that one is well in the hydrodynamic regime at intermediate temperature. In Fig. [21 we plot the the temperature 
dependence of the relative weights for g' = 0.5 and g' = 1.02. We see that there is no significant difference between the 
two results. In both cases, the first sound mode is dominant at low temperature, while the second sound is dominant 
at high temperature. This can be understood as follows. In the case of the supcrfluid velocity- velocity correlation, the 
external perturbation is directly coupled to the condensate motion. Therefore, the condensate mode should be always 
dominant in the spectral weight in the case of weakly- interacting Bose gas. At very low temperature (near T = 0), 
the first sound mode is essentially the condensate collective mode and the second sound mode is the collective mode 
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of quasiparticlc excitations. With increasing temperature, there is a crossover between two modes, and the nature of 
the sound oscillations changes. At high temperatures, the first sound mostly involves the noncondensatc oscillation, 
while the second sound mostly involves the condensate oscillation. In Appendix A, we compare the result with the 
self-consistent Hartree-Fock(HF) approximation, we see that the qualitative behaviors are well captured by the HF 
approximation. Moreover, we can explicitly see that in a weakly-interacting Bose gas, the single-particle spectrum is 
dominated by the condensate mode. 
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FIG. 2: Temperature dependence of relative weights of first and second sound modes Xi and X2 for g' — 0.5 and g' = 1.02. 

In Fig. [3j we compare the relative weights of first and second sound modes in both the single-particle spectral 
density A(q, lo) and the dynamic structure factor 5(q, uj), where 

q f Z 



where 



and 
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(48) 
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(50) 



We see that in S((\,uj), the hrst sound mode is dominant at all temperatures. This is in sharp contrast with A(q,uj), 
where the second sound mode is dominant near T c . 

IV. SINGLE-PARTICLE SPECTRAL DENSITY 



We now discuss the single-particle spectral density in connection with the experiment such as the the photoemission 
spectroscopy jl^ . In a uniform gas, the photoemission current is related to spectral weight 7 sw (q, w) cx A(q, w)/(oj) 
[lol [l^ . In Fig. 21 we show |A(q, w)/(o;)| for g' = 0.5 at the temperature T/T° = 0.5 as a function of q and w. In 
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FIG. 3: Comparison of the relative weights of first and second sound modes between A(q,u)) and S(q, co). 



the plot, the energy delta function is replaced by a gaussian with a finite width: 5(ui — c\q) — > Ae a ^~ Ciq ^ , where 
a = l/4dq 2 and A — \J~^-- Some appropriate value of the width parameter d is used. In Fig. 21 the wavenumber q 
and frequency uj are normalized in terms of the condensate healing lenght £ and the mean- field frequency ujq = gn/h. 
For a parameter set we used, £ « 0.87/im and ujq « 2.7kHz. The wavenumber and frequency are comparable to those 
in the recent experiment 
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FIG. 4: Plot of A(q,uj) as a function of q and u) for g' = 0.5 T = 0.5T C °. 



For comparison, we also calculate ^4(q, w) in the collisionless limit using the HFB-Popov approximation: 

A p °P° v (q,u) = u 2 q 5{u - E q ) - v 2 q 5{u + E q ) (51) 

In Fig. [5J we plot this collisionless result for the same coupling and temperature g' ~ 0.5, T — 0.5T C °. In the 
collisionless limit, there is only one sound mode. Moreover, it is clear from the expression (|5"Tj) that the weights of 
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uj = +E q and uj = —E q arc different, since in general u q ^ v q . This is in contrast with the two-fluid hydrodynamic 
regime, where uj = +1*1,29 and uj = —1*1,23 have the same weights. In particular, for large q one has u q 3> v q and thus 
A(q,uj) only has the contribution from uj = +E q . In the opposite low-frequency limit, one has u^/Vp — > 1. 

\A(q,co)f m (co)\ 
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FIG. 5: Plot of A Popov (ci,Lj) in the collisionless limit for g' = 0.5, T = 0.5T C . 



V. CONCLUSION 

In this paper, we have discussed single-particle spectral densities of first and second sound in supcrfluid Bosc gases. 
In order to obtain the thermodynamic quantities available for calculating the single-particle spectral densities of the 
first and second sound modes, we use the HFB-Popov approximation. We showed that that both first and second 
sound mode can be observed by photoemission spectroscopy near T ~ 0.5T C . We showed that the first sound mode is 
dominant at low temperature, while the second sound is dominant at high temperature. With increasing temperature, 
there is a crossover between two modes, and the nature of the sound oscillations changes. We hope that our results will 
stimulate further experiment by photoemission spectroscopy in a supcrfluid Bose gas in the two-fluid hydrodynamic 
regime. 

We have compared the relative weights of first and second sound modes in both the single-particle spectral density 
A(q, uj) and the dynamic structure factor S(q, uj). We showed that in 5(q, uj), the first sound mode is dominant at 
all temperatures. This is in sharp contrast with A(q, oj), where the second sound mode is dominant near T c . 

For illustration, we also considered the self-consistent Hartrcc-Fock approximation. We showed that the qualitative 
behaviors are well captured by the HF approximation, but the quantitative details are quite different. 

Finally, we showed thespectral weight \A(q,uj)f(uj)\ as a function of q and uj. We found that both first and second 
sound mode can be observed by photoemission spectroscopy. 
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Appendix A: Self-consistent Hartree-Fock approximation and the ZNG hydrodynamics 

For illustration, we consider the self-consistent Hartree-Fock approximation. In this approximation, the noncon- 
dcnsatc density is given by 

n = j^93/2(z), (Al) 



where A is the thermal de Broglie wavelength defined by 



2rf \ 1/2 



A = Us^J ■ < A2 > 

and the fugacity z is given by 

z = exp(-/3gn c ), (A3) 

where n c = n — n. The pressure is given by 

P = P + 7^g(n 2 + Inn - h 2 ) = P + -xg{n 2 c + 4nn c + 2n 2 ), (A4) 



where P is the kinetic pressure given by 

P — 

A 3 



P=^95/2(z) (A5) 



S = l(^P + gnn c ) (A6) 



The entropy is given by 

.<? = 

T \2 

We can then calculate the various thermodynamic derivatives used in the single-particle spectral function. 

Alternatively, one can use the linearized ZNG hydrodynamic equations [l|, [l8[ , which are based on the self-consistent 
HF approximation, to derive the supcrfluid velocity-velocity correlation function. The coupled equations for the 
condensate and noncondensate variables are given by 

^ = -n (V-5v„) + £r 12 , (A7) 
at 

mh —^- = -V5P -2gh Q \7(8n + Sn c ), (A8) 

dSP 5 ~ , r N - 2 , x , , . s 

= - T} Po(V ■5v n )-VP a + -{u cX) -U )5r 12 , (A9) 

-n c0 (V • £v c ) - ST 12 , (A10) 



dt 3 uv n/ u 3 
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dt 

m ^r = - vs ^ ( An ) 

at 

where 

<5/i c (r, t) = gSn c (r, t) + 2g5n(r, t) . (A12) 

The expression for <5r 12 is given by 

£Ti 2 [/] = - P° ncQ SfidiS, Mdiff = A - Mc- (A13) 

Tl2 
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In order to calculate the supcrfluid velocity-velocity correlation function, we add a time- dependent external current 
<5jex(r,i) to the equation of motion for 5n c given in (|A10|) . That is, we use 

d8n r 



at 



= -n c0 (V • 5v c ) - V • <5j ex - <5ri2, 



(A14) 



To solve the hydrodynamic equations, we introduce velocity potentials according to 5~v c = V</> c and <5v„ = V</>„. 
In terms of these new variables, the equations for the condensate and the equations for the noncondensate can be 
combined to give 

d 2 (t>c 2 2 a H 

gn c0 V 0c + 2gn V (f> n H o/idiff + ffV • oj cx , 
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Here <5ri2 has been expressed in terms of S^tdiS using (|A13[) . The equation of motion for S^dis is given by (15.71) 

(A17) 



dSfAdiff 2 2 2 5/idiff 
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dt 

We now consider an external current which excites modes of frequency u> and wavevector q 

5jex(r,t)=5jex(q,w)e i(tI - r - wt) , 
and look for the plane-wave solutions $ cn (r,t) = (/>c,n,<i,u^' i ' T ^ u3t - In this case, (|A17|) reduces to 



SfJLdiS 



1 — iuiT u 



gn C Q 



- ^<f> n ) q 2 -igq-S} e 



(A18) 



(A19) 



Substituting this result into (|A15I) and (|A16[) . we are left with two coupled equations for the supcrfluid and normal 
fluid velocity potentials: 



mil) <p c , 



gn-co 1 - 



-iff 1- 



<TH 



1 - iojT, 



1 - iuiT,, 



I q 2 <t>c^ + 2gh 



|q-5jex(q,w), 



oh n c0 
3(1 - iujt^) n 



(A20) 



and 



— 3- + 2pn 
3n 



2gn c0 



1 - 



2gg rtg 
9(1 - iwr^) rip 



i2g 



3(1 - iuJT u ) n Q 
oh n c0 



3(1 - iwr^) n 



q • 5j ex (q,w). 



(A21) 



Taking the limit ujt^ — > of these coupled equations, we obtain 



mu} 2 <f), 



ffn c0 (l - cr_ff)g 2 c , q ,w + 2.gn ( 1 + 



OHn c o\ 2 



-ig(l - o H )q • (5j ox (q,w), 



raw 



2ggng 
9ng 



q 4>r, 



+ 2gn c0 1 + 



ffgWeO \ „2 

3n 



9 ^c,q,w - «2.g 1 



OHn c o 
3n 



q- <5j cx (q,w). 



(A22) 



(A23) 



13 




14 



where 



Xi = 



x 2 = 



u 2 (uj - v\) 



(A35) 



The above expressions should be compared with the general expressions If we use the expression for u 2 to order 
g 2 given by (15.86), one has 

2 _ gn c0 imfio 4#n 



v 2 =u 2 + Av z , Av 2 
and the weights X\ and X 2 can be approximated as 

^-1 ~ — n n, Al 



"i 5P 



(A36) 



1 - 



"2 



1,^ Lit) 

,2 _ „,2„2 



(A37) 

u^q 2 is larger than the contribution at the 



"1 "2 

We find from (|A37[) that the contribution at the second sound mode w 
first sound mode u 2 = u\q 2 . 

Figure |6] compares the sound velocities calculated by the HF approximation with those calculated by the HFB- 
Popov approximation. Figure [7] compares X\ and X 2 calculated by the HF approximation with those calculated by 
the HFB-Popov approximation. We see that the qualitative behaviors are well captured by the HF approximation, 
but the quantitative details arc different. 




FIG. 6: Sound velocities of the first sound and second sound modes, dash Lines show the results from the self-consistent HF 
approximation, solid lines are the results from the HFB-Popov approximation. g' = 0.5. 



Appendix B: Validity of Hydrodynamics 

Here we discuss the validity of using Landau's two-fluid hydrodynamics in the experiment of Ref. |6j. In Fig. [51 
we plot the temperature dependence of the collisional relaxation time r M defined in Ref. [Lsl . [l|J] . This relaxation 
time describes the rate of equilibration of the condensate and noncondensate chemical potential. The condition for 
the hydrodynamic regime is described by ujt^ < 1 where uj is the frequency of the collective mode. We see that the 
collective mode of the frequency u ~ gn/h is well within the hydrodynamic regime in the intermediate temperature 
region T ~ 0.5T C . 



15 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



T/T 

c 

FIG. 7: The temperature dependence of the amplitudes Xx,X2- dash Lines show the results from the self-consistent HF 
approximation. Solid lines are results from the HFB-Popov approximation, g — 0.5. 
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FIG. 8: Temperature dependence of the relaxation time r^. g' = 0.5. 
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